Unified nonequifibrium dynamical theory for exchange bias and training eff^ects 
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We investigate the exchange bias and training effects in the FM/AF heterostructures using a 
unified Monte Carlo dynamical approach. This real dynamical method has been proved reliable and 
effective in simulating dynamical magnetization of nanoscale magnetic systems. The magnetization 
of the uncompensated AF layer is still open after the first field cycling is finished. Our simulated 
results show obvious shift of hysteresis loops (exchange bias) and cycling dependence of exchange 
bias (training effect) when the temperature is below 45 K. The exchange bias fields decrease with 
decreasing the cooling rate or increasing the temperature and the number of the field cycling. With 
the simulations, we show the exchange bias can be manipulated by controlling the cooling rate, the 
distributive width of the anisotropy energy, or the magnetic coupling constants. Essentially, these 
two effects can be explained on the basis of the microscopical coexistence of both reversible and 
irreversible moment reversals of the AF domains. Our simulated results are useful to really under- 
stand the magnetization dynamics of such magnetic heterostructures. This unified nonequifibrium 
dynamical method should be applicable to other exchange bias systems. 

PACS numbers: 75.75.+a.75.20.-g,75.60.-d,05.70.Ln 



INTRODUCTION 



Usually, when the heterostructure consisting of cou- 
pled ferromagnetic (FM) and antiferromagnetic (AF) lay- 
ers is cooled in field below the Neel temperature of its 
AF component, it shows the asymmetric magnetization 
P, 0, S 0, H, 01, which is referred to as the exchange 
bias effect. Furthermore, the exchange bias field, defined 
as the average of the two coercive fields, is observed to 
decrease with increasing the number of the consecutive 
field cycling, which is referred to as the training effect 
The exchange bias and training effects are very interest- 
ing and could be used in future spintronicsfsl. [a. [lOj and 
data storage. Usually, the FM layer is taken as a whole 
and the AF layer consists of many grains. The AF grain 
is small enough to consists of a single domain, and some 
uncompensated domains (or grains) may be formed by 
defects or impurities [ill, [l^ . [isj and couple with each 
other and with the FM domains. As the heterostruc- 
ture is cooled to a low temperature, the uncompensated 
spins in the grains and domains become locked-in and 
prefer to a unidirection in the interface, thus contribute 
to the magnetization shift Moreover, under the rever- 
sal of FM domains, the uncompensated grains or domains 
will be irreversibly reorganized [13, [l^ 0, [ill and 
thus cause the training effect. The idea of domain states 
was corroborated in some Monte Carlo simulations [Toll. 
On the other hand, Hoffmann 20] considered the biaxial 
anisotropy of the AF sublattices and solved it by varia- 
tional method. Actual nonequifibrium dynamical prop- 
erties of the magnetization are still waiting to be eluci- 
dated. It is highly desirable and needed to systematically 
investigate the two effects in a unified theory. 

In this article we use a unified Monte Carlo dynam- 
ical approach [U to study the FM/AF heterostructure 
in order to investigate the exchange bias and training 



effect. Our simulated result shows the obvious shift of 
hysteresis loops and the cycling dependence of exchange 
bias. The magnetization of uncompensated AF layer is 
still open after the field cycling is finished. The exchange 
bias fields decrease with decreasing the cooling rate or 
increasing the temperature and the number of the field 
cycling. With the simulations, we shows the exchange 
bias can be manipulated by controlling the cooling rate, 
the distributive width of the anisotropy energy, or the 
magnetic coupling constants. Essentially, these two ef- 
fects can be explained on the basis of the microscopically 
irreversible reversal of the AF domains. More detailed 
results will be presented in the following. 

The remaining part of this paper is organized as fol- 
lows. In next section we shall define our model and 
discuss our simulation method. In section III we shall 
present our simulated results and analysis. In section 
IV we shall discuss the microscopic mechanism for the 
phenomena in a unified way. Finally, we shall give our 
conclusion in section V. 



II. MODEL AND METHOD 

According to experimental observations [22j . for both 
compensated and uncompensated AF layers the easy axis 
tends to form along external cooling field direction rather 
than later rotating field direction. In our model the AF 
layer consists of many AF domains, and the FM layer 
consists of one single domain. Assuming the cooling field 
is applied parallel to the AF/FM interface, then all the 
easy axes of AF and FM domains lie in the plane of the 
interface. The coupled bilayers of AF and FM domains 
are shown in the inset of Fig. 1(a). The rectangles of the 
white pattern represent the AF domains and the larger 
rectangle is the single FM domain. The AF domains cou- 
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pie to each other antiferromagnetically and the single FM 
domain couples to all the AF domains ferromagnetically. 
We define the z axis along the common easy axis which 
lies in the interface plane. We apply the external field to 
saturate the magnetization of the FM layer along the z 
axis. 

For simplicity, we consider all the uncompensated spins 
in the AF domains are the same. We use S'si to denote 
the spin vector of the jth AF domain and Ss to denote 
that of the single FM domain, where S' and S are the un- 
compensated spin values and FM spin respectively. Then 
we write the Hamiltonian of the bilayers in an external 
field as 

i i 

+ ^1 ^ Si • - J2 ^ ■ s (1) 
id i 

where 7' = g^ia^sS' and 7 = qiiq^bS. The first and sec- 
ond terms represents the anisotropy of the FM domain 
and the AF ones, and Ku and kui are the corresponding 
anisotropy constants. The third term represents the Zee- 
man energy of the moments due to the applied external 
field. The fourth term represents the antiferromagnetic 
coupling among the AF domains. The last term repre- 
sents the ferromagnetic coupling between the FM and AF 
domains. 

Using 9i and (3 to describe the angles of the i-th. AF 
moment and the FM moment deviating from the common 
easy axis, we can express the energies of the FM domain 
and the i-th AF as 

H^^ = ~{J2^SiS + Kucos 13 + -fBs) COB (3 (2) 

i 

and 

Hf'^ = {JiSi'^Sj-J2SiS-kuiCOs9i--i'Bsi)cos9i (3) 

where both Si and s are the scalars taking either 1 
or -1. Thus for the i-th AF domain the energy in- 
crement is Ai?i = kuisiv? 9i — hi{cos9i — 1), where 
hi = (— Ji Sj + J2S + j'B)si, and for the FM domain 

the energy increment is AE = Ku sin^ /3 — hp (cos (3 — 1), 
where hp = (J2 J^i + lB)s. We can express IS.E and 
AE, asil 



and 
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As a result, to reverse its moment, the the FM layer 
must overcomes a barrier E^ = K^i^ + hp/2Ku)'^ if 
\hp\ < 2Ku, or 2hp if hp > 2Ku\ and the i-th AF grain 



a barrier El = kui{l + hi/2kui)'^ if \hi\ < 2fc„i, or 2hi if 
hi > 2kui. If the condition hp < —2Ku or hi < — 2fc„,j is 
satisfied, there is no barrier for the reversal. 

Actually, for the distribution of the AF anisotropy en- 
ergy we use a Gauss function, f{kui) — exp[— (fc„i — 
kuY /cr"^], whose a and ku are set to 30.0 meV and 50.0 
meV unless stated otherwise. The anisotropy energy of 
the FM domain is set 200.0 meV without losing main 
physics. Thus the reversal rate for a spin to reverse is 
R = RQe~^''/^'^^ , where Ei, is the energy barrier and 
i?o is the characteristic frequency. In our simulations, 
i?o is set to 1.0 X 10^/s. We adopt a square lattice for 
the AF domains and use 20 x 20 as its size. Since we are 
only interested in the exchange bias and training effect at 
the nanoscale, the AF lattice is enough to capture main 
physics. Furthermore, we assume the AF domains have 
uniform moment 4.0 /is and the FM domain 2000 /is. 
The coupling constant Ji is set to 4.0 meV, and J2 8.0 
meV. In our simulations the system is quenched from a 
high-enough temperature such as 610 K, at which the AF 
layer is paramagnetic, to a low-enough temperature such 
as 10 K. The magnetization and exchange bias fields are 
calculated at the low temperature 10 K unless the tem- 
perature is explicitly stated otherwise. The basic rate of 
changing temperature is vq= 50 K/s. The field sweeping 
rate is set to 0.5 T/s with the basic increment 0.1 T for 
each simulation step. 



III. SIMULATED RESULTS AND ANALYSIS 

At first, we let the AF/FM bilayers relax under a mag- 
netic field of 5.0 T at a high temperature 610 K. This 
temperature is enough to make both the FM layer and 
the AF layer remain paramagnetic. When the tempera- 
ture decreases, the average magnetization values of the 
two layers increases. The external field makes the aver- 
age magnetization of the FM layer have a large increase 
below 600 K, and reach nearly to the saturated value at 
500 K. When the temperature becomes lower than 60 
K, the average magnetization of the AF layer looks like 
that of an antiferromagnet under an applied field and is 
dependent on the cooling rate v. Then, we further cool 
the bilayers under the same field. After the temperature 
reaches down to 10 K, we start to change the field while 
keeping the temperature unchanged. The field decreases 
from 5.0 T to -10.0 T and then increases back to 5.0 T for 
the first hysteresis. Repeating the field cycling, we will 
make the second hysteresis loop. The simulated results 
are shown in Fig. 1. 

As shown in Fig. 1(a), the origin of the first hysteresis 
is clearly shifted in the negative field direction and shows 
the exchange bias. The exchange bias field is defined as 
Hp = {Hci + Hcr)/2, where Hd and Her is the coercivity 
of the left and right branches. The left branch of the sec- 
ond hysteresis moves towards the positive direction, but 
the right branches of the first two loops almost coincide 
with each other. Actually, any further loop almost does 
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FIG. 1: The first two hysteresis loops of the FM (a) and AF 
(b) layer at 10 K. The inset in (a) shows the AF/FM bilayers. 
The hysteresis loop is obtained by changing the field in the 
order of P1-P2-P3-P4-P1. 



no difference in the right branch from the second hys- 
teresis. The shift of the second loop clearly demonstrates 
that the bilayers magnetization depends on the cycling 
history, which is known as training effect. Fig. 1(b) 
shows the magnetization of the AF layer, which drops 
largely and opens widely due to the irreversible reversal 
of the AF domains after the first field cycling is finished. 
The subsequent magnetization is more smooth but still 
not closed, indicating the continuing cycling dependence 
of exchange bias. This is consistent with other Monte 
Carlo simulations . 

We study the effect of the temperature T on the ex- 
change bias field. He, for different loops. Our simulated 
exchange bias fields as functions of T for the first two 
loops are shown in Fig. 2. For both of the two curves, 
the data can be fitted by the simple function 



(6) 



where ai, Ti, and bi are fitting parameters. For the 
fitting in Fig. 2, the parameters ai, Ti, and bi takes 
4.26 T, 17.95 K, and 1.58 for the first loop, and 3.94 
T, 16.66 K, and 1.59 for the second loop. Our results 
are consistent with experimental observation that the ex- 
chaiige bias field decreases with increasing temperature 

The exchange bias field is dependent on the field cy- 
cling number n. Our simulated result from n=l to n=9 
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FIG. 2: Temperature dependence of the exchange bias fields 
for the first two loops. The exchange bias field is calculated 
at a given temperature after the system is cooled from 610 K 
to the temperature value. The lines are the fitting curves in 
terms of the simple function defined in Eq. (6). 




FIG. 3: The loop-number dependence of the exchange bias 
fields for the Gaussian width (t=30.0 and 20.0 meV. The tem- 
perature is 10 K. All the data except for n — 1 can be well 
fitted by a simple function — ^qH E(n) — a2p" + 62. 



is shown in Fig. 3. Here, the temperature is kept at 10 K, 
and (7 is set to 20.0 and 30.0 meV. For both of the curves 
in Fig. 3, the data points excepts of n = 1 are well fitted 
by the simple function —fioHEiji) = a2p" -I- 62, where 
02, 621 and p are the fitting parameters, taking 0.23 T, 
2.26 T, and 0.76 for a = 30.0 meV, and 0.22 T, 1.51 T, 
and 0.74 for a = 20.0 meV. The value of -^moHe{1) usu- 
ally is substantially above the extrapolation of the other 
—lii)HE{n) [n > 1). These simulated results are in good 
agreement with experimental observation [7]. 

Since the training effect reflects the non-equilibrium 
dynamical magnetization which is caused by the irre- 
versible reversal of meta-stable domains formed during 
quenching, the quenching rate must affect the exchange 
bias field. By changing the cooling rate we study the 



4 




FIG. 4: The cooling-rate dependence of the exchange bias 
fields for the first two loops. The lines are the fitting curves 
in terms of Eq. (7) . 



exchange bias field as the function of quenching rate v. 
The result is shown in Fig. 4. For both of the loops, the 
data be well fitted by 



IMqHe = as \n{b3 h 1) 

1^0 



(7) 



where and 63 are 0.292 T and 20847 for the first loop, 
and 0.262 T and 16248 for the second loop. The ex- 
change bias field at 10 K increases logarithmically with 
increasing the quenching rate. 

It is interesting to investigate the dependence of the 
exchange bias field on the coupling constants Ji and J2. 
The simulated results are shown in Fig. 5. As shown in 
Fig. 5(a), the exchange bias field decreases with increas- 
ing Ji . The training effect is nearly unchanged when Ji 
changes from to 2meV, but diminishes to zero quickly 
with increasing Ji from 2meV. When Ji is larger than 6 
meV, the exchange bias field already becomes very small 
and the training effect is actually zero. In contrast, the 
J2 data points in Fig. 5(b) can be well fitted by the 
simple function — /iq-H^e = a4(exp( J2/64) — 1), where 04 
and 64 are the fitting parameters, taking 0.44 T and 4.01 
meV for the first loop, and 0.51 T and 4.49 meV for 
the second loop. It is clearly shown that the exchange 
bias field increases exponentially as J2 increases. Both 
of the the exchange bias field and the training effect can 
be enhanced by decreasing Ji and increasing J2, or by 
increasing J2/J1. This is consistent with experimental 
trend [2^. In Fig. 6 we shows how the distributive width 
a of the AF anisotropy affects the exchange bias fields. 
Clearly the exchange bias field increases with a, and so 
does the training effect. Experimentally, the width can 
be increased by the additional nonmagnetic impurities 
and the enhanced roughness of the AF crystalline phases. 
This implies that the rougher the AF crystalline phases 
are, the larger the exchange bias and training effect. Our 
result reveals that the exchange bias field is determined 
by both the coupling constants and the distributive width 
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FIG. 5: The exchange bias fields as functions of the coupling 
constants Ji (a) and J2 (b) for the first two loops. 
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FIG. 6: The exchange bias field as a function of the Gaussian 
width a for the first two loops. 



of the AF domain anisotropy. These are useful to com- 
pletely understand the phenomenon j25l. [26j . 



IV. TRENDS AND MICROSCOPIC 
MECHANISM 

After being cooled down to the low temperature, the 
AF layer has a non-zero net FM moment due to the 
driving of both the field and the FM layer. Assuming 
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there are Na AF domains, on average we have the mo- 
ments in part of all the AF domains aligning parallel 
although they are coupled with AF interactions. The ex- 
change bias field is determined by the effective moment 
Ma, the difference of Ma between the first two loops 
determines the training effect. Naturally, both the ex- 
change bias field and the training effect increase with 
increasing J2 and with decreasing Ji, as shown in Fig. 5. 
Actually, small Ji does not affect the effects, but larger Ji 
than 2meV is harmful to the effects at 10 K. In addition, 
it is easily understood that Ma decreases with increasing 
the temperature T. As a result, both the exchange bias 
field and the training effect decrease with increasing T, 
as shown in Fig. 2. The exponential description in Eq. 
(6) refiects the fact that moment reversals are thermally 
activated. It is reasonable that both the exchange bias 
field and the training effect increase with increasing the 
cooling rate i', as shown in Fig. 4. This is mainly because 
the average magnetization of the AF layer increases with 
v when the temperature is below 50 K. When the cooling 
rate approaches to zero, both the exchange bias field and 
the training effect should be zero. In another word, our 
results should approach to those of corresponding equi- 
librium systems when the cooling rate v approaches to 
zero. 

As shown in Fig. 6, both the exchange bias field and 
the training effect are nearly zero when the distributive 
width a of the AF anisotropy energy is smaller than 15 
meV, but they increase substantially with increasing a 
for (T > 15 mcV. This means that the effects are depen- 
dent on a wide distribution of the anisotropy energy. This 
can be understood in terms of the changing of the energy 
barriers with the external field. From PI to P2 in Fig. 
1(a), the effective barrier of the FM layer decreases but 
is still high enough to avoid the reversal, but meanwhile, 
more and more spins of the AF domains are reversed due 
to their lower energy barriers. At the point P2, the FM 
moment is reversed with the help of the field and the re- 
versing of the AF domains. Anyway, some AF domains 
with high energy barriers have their moments unchanged, 
even after the FM layer has been reversed, and thus there 
is a net average moment of the AF domains parallel to 
the moment of the FM layer. This net average moment 
increases with the distributive width a. This explains the 



increasing of the exchange bias field and training effect 
with increasing a. The more the field cycling loops, the 
longer the time. Actually, this is similar to reducing the 
cooling rate in effect. As a result, the exchange bias 
field decreases with increasing the number of the field 
cycling. The turning point of the time scale causes the 
largest drop happens between the first loops. 



V. CONCLUSION 

In summary, we use a unified Monte Carlo dynami- 
cal approach |2ll| to study the FM/AF heterostructure in 
order to investigate the exchange bias and training ef- 
fect. The magnetization of uncompensated AF layer is 
still open after the first field cycling is finished. Our sim- 
ulated result shows the obvious shift of hysteresis loops 
(exchange bias) and the cycling dependence of exchange 
bias (training effect). The exchange bias fields decrease 
with decreasing the cooling rate or increasing the tem- 
perature and the number of the field cycling. With the 
simulations, we show the exchange bias can be manip- 
ulated by controlling the cooling rate, the distributive 
width of the anisotropy energy, or the magnetic coupling 
constants. Essentially, these two effects can be explained 
on the basis of the microscopical coexistence of both re- 
versible and irreversible moment reversals of the AF do- 
mains. Our simulated results are useful to really under- 
stand the magnetization dynamics of such magnetic het- 
erostructures which should be important for spintronic 
device and magnetic recording media [25l . [26l . [27j . This 
unified nonequilibrium dynamical method should be ap- 
plicable to other exchange bias systems. 
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